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ABSTRACT 

We present photometry of WASP- 10 during a transit of its short-period Jovian planet. We employed 
the novel PSF-shaping capabilities of the OPTIC camera mounted on the UH 2.2 m telescope to 
achieve a photometric precision of 4.7 x 10~ 4 per 1.3 min sample. With this new light curve, in 
conjunction with stellar evolutionary models, we improve on existing measurements of the planetary, 
stellar and orbital parameters. We find a stellar radius R* — 0.698 ± 0.012 i? Q and a planetary 
radius Rp = 1.080 ± 0.020 Rj up . The quoted errors do not include any possible systematic errors in 
the stellar evolutionary models. Our measurement improves the precision of the planet's radius by a 
factor of 4, and revises the previous estimate downward by 16% (2.5a, where a is the quadrature sum 
of the respective confidence limits). Our measured radius of WASP- 10b is consistent with previously 
published theoretical radii for irradiated Jovian planets. 

Subject headings: instrumentation: detectors — techniques: photometric — stars: individual (WASP- 
10) — planetary systems: individual (WASP-lOb) 
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1. INTRODUCTION 

When an exoplanct transits its star it provides a 
valuable opportunity to study the physical character- 
istics of planets outside our solar system. Photo- 
metric survey teams are discovering a rapidly grow- 
ing number of short-period planets transiting relatively 
bright (9 < V < 13) Sun-like stars by monitoring 
the brightnesses of hundreds of thousands of stars for 
weeks or months, and following up on transit can- 
didates with spectroscopy to rule out false positives 
(lAlonso et alJl2004t iMcCullough eta l. 2006: Bakos et all 
120071 ICameron et all l2007t iBarge etall I2008D . With 
this windfall of new planet discoveries comes a need 
for efficient and precise photometric follow-up. Mea- 
suring planet properties precisely enough for meaning- 
ful tests of theoretical models of planetary interiors re- 
quires transit light curves with good photometric pre- 
cision (<10~ 3 ) and high cadence (<1 min). In the 
past these requirements have often been met by combin- 
ing photometry from ma ny transit events (jHolman et al.l 
2006J IWinn et al.ll2007bl ). or observing from space using 
the Hubble and Spitzer space telescopes (|Brown et al.1 
120011: iGillon et all 120071) . These approaches work but 
they are costly in time and/or resources. 

Our approach is to capitalize on the strengths of a rela- 
tively n ew type of detector : the orthogonal transfer array 
(OTA; iTonrv et all 1 19971) . OTA detectors are charge- 
coupled devices (CCDs) employing a novel four-phase 
gate structure to allow accumulated charge to be shifted 
in two dimensions during an exposure. Usually, this ca- 
pability is used to compensate for image motion, pro- 
viding on-chip tip-tilt correction and a narrower point 
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spread funct i on (P SF). However, as demonstrated by 
IHowell et all (|2003f) . OTAs can also be utilized to delib- 
erately broaden PSFs in a manner that is well suited for 
high-precision photometry. The method of "PSF shap- 
ing" uses a raster scan to distribute the accumulated sig- 
nal from each star over many pixels within a well-defined 
region of the detector, thereby collecting more light per 
exposure and increasing the observing duty cycle com- 
pared to traditional CCDs. This OTA observing mode 
provides benefits similar to defocusing, but without the 
downside of a spat i ally and chromati cally varying PSF 
(IHowell et al.ll2003 ITonrv et al.ll2005l ). 

In this Letter, we present new OTA-based photometry 
for a transiting planetary system recently announced by 
the Super Wide- Angle Search for Planets (SuperWASP) 
consortium. The planet, WASP-lOb, is a 3.3 Mj up planet 
in an eccent ric, 3.09 day orbit ar ound a V = 12.7 Re- 
type dwarf (jChristian et al.l 120081 . hereafter C08). C08 
measured a radius for WASP-lOb that is much larger 
than predicted by models of planetary interiors. In §[2]we 
present our observations. We describe our data reduction 
procedures in §[3] and our light curve modeling procedure 
and error analysis in § HJ We present our results in § [5] 
and conclude in §[6]with a brief discussion of our findings. 

2. OBSERVATIONS 

We observed WASP-10 on UT August 16, spanning a 
transit time predicted by the ephemeris of C08. We used 
the Orthogonal Parallel Transfer Imaging Camera (OP- 
TIC) on the University of Hawaii 2.2 m telescope atop 
Mauna Kea, Hawaii. OPTIC consists of two 2048 x 4096 
pixel Lincoln Lab CCID28 OTA detectors with 15 /xm 
pixels, corres ponding to 0*7135 p ixel -1 and a 9'2 x 9'2 
field of view (|Howell et al J 12003). Each detector is read 
out through two amplifiers, dividing the frame into four 
quadrants. The portions of each quadrant closest to the 
amplifiers can be read out without disturbing the rest of 
the detector. Small subregions near the amplifiers can be 
read out rapidly (at a rate of up to 100 Hz) to monitor 
the position of one to four guide stars. The centroids 
of the guide stars are analyzed and can be used to shift 
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Fig. 1. — Example of an OPTIC PSF selected from the WASP-10 field. Throughout the exposure, the accumulated signal from the 
incident stellar PSF is shifted into neighboring pixels to produce a PSF confined to an approximately box-shaped region. The inset in the 
lower left region shows a separate, 5-second "snap-shot" exposure of WASP-10 with no OT guiding. The solid line shows the star aperture 
used to calculate the instrumental magnitude of the star. The background level is estimated from the mean level measured from the four 
flanking boxes (dashed lines). 



charge in the science regions to compensate for image 
motion due to seeing and telescope guiding errors. In 
our application, we used a single OTA guide star as a 
reference point to define the origin of the shaped-PSF 
raster scan. We maximized our field of view by using the 
remaining three guide regions to image the target field. 

We observed WASP-10 continuously for 4.35 hr span- 
ning the predicted mid-transit time. Conditions were 
photometric. The seeing ranged from 0'.'5 to 1'.'2. We ob- 
served through a Sloan z 1 filter, the reddest broad-band 
filter available, to minimize the effects of differential at- 
mospheric extinction on the relative photometry, and to 
minimize the effects of limb darkening on the transit light 
curve. We shifted the charge during each integration to 
form a square PSF with a side length of 20 pixels (2'.'7). 
During each exposure, the on-chip raster scan was traced 
out 11 times to uniformly fill in the shaped-PSF region. 
An example of a WASP-10 PSF is shown in Figure [TJ 
The integration time was 50 s, the readout time was 25 s, 
and the reset and setup operations took an additional 4 s. 
The resulting cadence was 79 seconds, corresponding to 
a 63% duty cycle. 

3. DATA REDUCTION 

To calibrate the science images, we first subtracted the 
bias level corresponding to each of the four amplifiers. 
Visual examination of a median bias frame revealed that 
the bias offset is nearly constant across the detector rows. 
We opted to estimate the bias level of the individual 



science frames using the median of each 32-pixel-wide 
over-scan row. For flat-field calibration, we normalized 
individual dome flats by dividing each quadrant by the 
median value. We then applied a median filter to the 
stacked set of normalized flats to form a master flat for 
each night. Finally, we constructed a customized flat for 
each science frame by convolving the unshifted master 
flat with the OTA shift pattern (which is r ecorded in an 
auxil iary file at the end of each exposure) (|Howell et al.l 
l200l . 

We performed aperture photometry using custom pro- 
cedures written in IDL. The flux from each star was mea- 
sured by summing the counts within a square aperture 
centered on the PSF. Before the summation, we sub- 
tracted the background level, taken to be the (outlier- 
rejected) mean count level among four rectangular re- 
gions flanking the star aperture. The star aperture and 
sky regions are illustrated in Figure [1] 

Ideally, the PSF shaping procedure confines the accu- 
mulated charge for each star within a well-defined region 
of the detector that is common to all images. However, 
the telescope focus gradually drifted, resulting in varia- 
tions in the full width at half maximum (FWHM) of the 
shaped-PSFs ranging from the nominal 20 pixels up to 
24 pixels before we noticed this problem and refocused 
the telescope. We also noticed a faint halo surrounding 
the main PSF, extending to several pixels beyond the 
nominal PSF width. To ensure that our digital aperture 
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Fig. 2. — Relative photometry of WASP-10. The solid red line shows the best— fitting model light curve and the residuals (O-C) are 
plotted belo. 



captures all of the incident flux, we used a very wide 
aperture width of 50 pixels, chosen because it gave the 
smallest scatter in relative flux outside of the transit. 
The sky regions extend from 30 to 70 pixels from the 
center of the PSF. 

To correct for variations in sky transparency, we di- 
vided the the flux of WASP-10 by the sum of the fluxes 
of 5 comparison stars with apparent magnitudes com- 
parable to WASP-10. We used the root-mean-square 
(rms) scatter of the out-of-transit (OOT) points as an 
estimate for the individual measurement uncertainties. 
Our relative photometric measurements of WASP-10 are 
shown in Figure [2J and listed in Table Q] along with with 
the heliocentric Julian dates (HJD) of observation. 

4. MODELING 

We determined the stellar, planetary and orbital char- 
acteristics by (1) fitting the light curve with a model 
consisting of a star and planet in Keplerian orbit around 
their center of mass; (2) combining the results with 
stellar-evolutionary models to break the modeling degen- 
eracies inherent in the light-curve model. 

4.1. Photometric Model 

Our light curve fit ting procedure is similar to that of 
iHolman et ail (|2006l ) an dlWinn et all (l2007aft. and uses 



the analytic formulas of lMandel fc Agoil ""( 2002f ) to com- 
pute the integrated flux from the uneclipsed stellar sur- 
face as a function of the relative positions of the star and 
planet. The parameters were the stellar radius i?* and 
mass M*, the planetary radius Rp and mass Mp, the or- 
bital inclination i, period P, eccentricity e, argument of 
periastron u>, and mid-transit time T c . We fixed P at the 



value 3.0927616 d given by C08, which has negligible un- 
certainty for for our purpose. Due to the well-known fit- 
ting degeneracies M p cx M^ 3 and R p oc i?* oc M^ 3 (see, 
e.g., Winn 2008) we fixed M* at a fiducial value, thereby 
solving effectively for i?*/a and Rp/a rather than i?* 
and Rp. The planet mass Mp is included for complete- 
ness but the data are highly insensitive to this parame- 
ter. We measured uj and e by fitting a Keplerian model 
to the radial velocity (RV) measurements reported by 
C08 and fixed the parameters in our subsequent analysis 
(e = 0.051, uj — 153 deg). Because the eccentricity is so 
close to zero, we found that the actual values had very 
little effect on our results. There were four other free pa- 
rameters: two parameters describing a linear function of 
airmass used to correct for differential airmass extinction, 
and the two coefficients in the quadratic limb darkening 
law. We fitted the model to our photometric observations 
using the goodness-of-fit statistic 

fijobs) - /j(calc) y 



* 2 = E 



where /^(obs) and /i(calc) are the observed and com- 
puted fluxes, respectively, and o~i is the uncertainty of 
each flux measurement, set equal to 0.000475 for our 
measurements, for reasons explained below. 

To determine the best-fitting parameter values and 
their uncertainties, we used a Markov Chain Monte Carlo 
algorit hm with a Gibbs sampler (see e.g. I Winn et al.l 
2008b, and references therein). We created 10 chains of 
10 6 steps each, and combined them after verifying that 
the chains were converged and well-mixed. The recorded 
sequence of parameter combinations provides the pos- 
terior probability distribution of each parameter. We 
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record the 15.9, 50.0, and 84.1 percentile levels in the cu- 
mulative distribution. For each parameter we quote the 
50.0% level (the median) as the "best value," and use the 
84.1% and 15.9% levels to define the "one sigma" upper 
and lower confidence limits. 

After optimizing the model parameters, we calculated 
the residuals to assess the photometric noise, and in par- 
ticular the time-correlated noise that may arise from in- 
strumental, atmospheric and astrophysical sources — the 
so-called "red nois e" — using the approach outlined by 
IWinn et alj (|2008bl . see also Gillon et al. 2006). We aver- 
aged the residuals into M bins of N points and measured 
the standard deviation, erjv, of the binned residuals. If 
the noise were described by a Gaussian distribution then 
the rms of each grouping of residuals should decrease as 
iV- 1 / 2 [M/(M-l)] 1 / 2 ox, where (7i is the rms of the un- 
binned residuals. Correlated noise would usually cause 
the observed scatter to be larger than this expectation, 
by some factor (3. We computed (3 for binning intervals 
ranging from 5 to 20 minutes, finding [3 ~ 1. For this 
reason we set our measurement uncertainties equal to 
a x = 0.000475 (0.5 mmag). 

4.2. Stellar Evolutionary Models 

2/3 

To break the fitting degeneracies M p oc MJ and 

R p oc R+ oc M*^ 3 , we required consistency between the 
observed properties of the star, the stellar mean den- 
sity p+ that can be derived from the photometric param- 
eter a/R+ (Seager & Mallen-Ornelas 2003, Sozzetti et 
al. 2007, Torres et al. 2008), and the Yonsei-Yale (Y 2 ) 
stellar evolution models by Yi et al. (2001) and Demar- 
que et al. (2004). The inputs were T eff = 4675 ± 100 K 
and [Fe/H] = 0.03 ±0.20 from C08, and the stellar mean 
density p+ = 3.099 ± 0.088 g cm~ 3 derived from the 
light curve. We computed isochrones for the allowed 
range of metallicities, and for stellar ages ranging from 
0.1 to 14 Gyr. For each stellar property (mass, radius, 
and age), we took a weighted average of the points on 
each isochrone, in which the weights were proportional 
to exp(^x*/2) with 
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Here, the A quantities denote the deviations between 
the observed and calculated values at each point. The 
weights were further multiplied by a factor taking into ac- 
count the number density of stars along each isochrone, 
assuming a Salpeter mass function. It is important to 
keep in mind that this procedure assumes that any sys- 
tematic errors in the Y 2 isochrones are negligible. 

5. RESULTS 

Table [2] lists the parameters for the WASP-10 plane- 
tary system derived from the analysis of our new pho- 
tometry. Of particular interest are the the stellar and 
planetary radii, i?» = 0.698 ± 0.012 R and R P = 
1.080 ± 0.020 -Rjup, and the mid-transit time, which is 
measured with a precision of 7 s. For completeness, we 
also list the orbit parameters determined from our inde- 
pendent analysis of the RV measurements of C08. Fi- 
nally, we list several other quantities derived from our 
light curves that can be used for additional follow-up 
observations and modeling efforts. 



Our measured planetary, stellar and orbital parame- 
ters differ significantly from the values reported previ- 
ously by C08. For example, our measurement of the 
stellar radius is smaller by 1.9<7, where a is quadrature 
sum of the respective confidence limits, and our mea- 
sured orbital inclination is higher by 2.7<r. As for the 
planetary radius, we refined the C08 value downward by 
2.5<7, while improving the precision from rj 0.08 i?j up 
to 0.02 i?j U p. The disagreement between the two results 
for the planetary radius can be partly attributed to dif- 
ferent stellar radii, as well as different transit depths: 
(Rp/R*) 2 = 0.02525±£^28 from our analysis, com- 
pared to {R P /R*f = 0.029 ± 0.001 from C08. 

We do not know the reason for these discrepancies. 
C08 reported four light curves: one based on the phased 
photometry from the 20 cm SuperWASP survey tele- 
scope, and three higher precision light curves observed 
with larger 1 m class telescopes. The latter three data 
sets were incomplete, covering ingress or egress but not 
both for a given transit. In our experience, complete 
light curves such as ours are much preferred, as they al- 
low for more thorough checking of the noise properties 
and any systematic effects. C08 did not report any tests 
for red noise or other systematic effects, although such 
effects are visually evident for the Tenagara light curve 
and may be present in the other light curves as well. In 
addition, C08 did not report whether they allowed the 
limb-darkening parameters to vary, or held them fixed at 
values deemed appropriate; it has been shown by South- 
worth (2008) and others that holding them fixed can re- 
sult in underestimated parameter errors. 

6. SUMMARY AND DISCUSSION 

We have presented high-precision photometry of the 
K-type star WASP-10 during the transit of its close-in 
Jovian planet. Our analysis has improved the precision 
with which the planetary, orbital and stellar properties 
are known. From the photometry, we improved the un- 
certainty in the orbit inclination angle by a factor of 2.7 
and the transit depth by a factor of 3.8, and we measured 
the transit mid-time to 7 seconds. We also measured the 
stellar density (/?*) from the transit light curve, which 
we used together with stellar evolution models to mea- 
sure the stellar mass to a precision of 4.6%. We used 
the refined stellar mass to improve the precision in the 
planetary radius by a factor of 4. Notably, the 1.8% un- 
certainty in Rp is dominated by the uncertainty in the 
stellar mass, rather than the photometry. 

Our determination of the planetary radius represents 
a 16% downward revision of the value reported by 
C08. Those authors found that the radius of WASP- 
10b was much larger than predicted by even a coreless 
planet model. Following C08, we compared our revised 
radius of WASP- 10b to the theoretical predictions of 
iFortnev et alJ (|2007l ). who provide tabulated radii for 
planets with a range of semimajor axes, masses, ages 
and heavy-element core masses. For a hydrogen- helium 
planet with properties appropriate to WASP-10, Fort- 
ney et al. predict a radius of 1.13 i?j up for the case of a 
coreless planet with an age of 4.5 Gyr. Adding 25, 50, 
and 100 M ffi of heavy elements (nominally in the form 
of a solid core) decreases the radius to 1.11, 1.09 and 
1.06 i?j U p, respectively. The models are able to repro- 
duce our measured radius (Rp = 1.08 ± 0.02 Rj U p) by 
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including a heavy-element core with a mass ranging from 
35 M e to 95 M®. 

We also compared our WASP- 10b mass and ra- 
dius measurements to the theoretical models of 
iBodenheimer et all (|2003l ). who predicted the radii of gas 
giant planets as a function of mass, age and equilibrium 
temperature. Assuming a planet temperature T cq = 
TeffjR Ja) 1 / 2 = 1370 (Table 2), IBodenheimer et all 
(2003) predict R P = 1.11 i? Jup for a 4.5 Gyr, solar- 
composition planet. The predicted radius is essentially 
the same whether the composition is taken to be purely 
solar or if a 40 M® core of heavy elements is included. 
It is also fairly insensitive to the value of T eq . Thus our 
me asured radius is only 3% (1 .5a) smaller than predicted 
bv IBodenheimer et all (|2003l ). 

We conclude that WASP-lOb is not significantly "in- 
flated" beyond the expectations of standard models of 
gas giant planets. Instead, we find that WASP-lOb is 
among the densest cxoplanets so far discovered. Its 
mass and radius similar to another transiting planet, 
HD 17156b (Winn et al. 2008c). However, WASP-10 
and HD 17156 are very different host stars. Whereas 
HD 17156 is metal-rich ([Fe/H]= +0.24) and has a mass 
of 1.26 M Q , WASP-10 is a 0.75 M© K dwarf with a 
much lower atmospheric abundance ([M/H]= +0.03±0.2; 
C08). WASP-10 stands out as one of only nine known 
low-mass (M* < 0.8 M ) planet host stars within 200 pc, 
and its planet is among the most massive so far detected 
in this sample of low-mass stars 5 . 

We have found the OPTIC instrument to be capa- 
ble of delivering single-transit photometry that is precise 
enough for meaningful comparisons of observed plane- 

5 http://cxoplanets.org 



tary properties and theoretical models. Our WASP-10 
light curve has a precision of 4.7 x 10~ 4 per 1.3 min sam- 
ple, or 6.8 x 10~ 4 per 1 min bin, for a V = 12.7 star 
(Figure [2]). This is comparable to the precision that has 
been reported for composite light curves, based on obser- 
vations of man y transits, with trad itional CCD cameras. 
For example, Uohnson et al.1 (|2008f ) presented a compos- 
ite z-band light curve for HAT-P-1 (V = 10.4) based on 
observations of 8 transits, with a precision of 5.7 x 10~ 4 in 
1 min bins. OTAs may also prove useful for the future de- 
tection and characterization of transiting Neptune-mass 
and "super-Earth" planets. 
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TABLE 1 

Relative Photometry for WASP-10 
Heliocentric Julian Date Relative Flux 



2454663.94419 
2454663.94512 
2454663.94604 
2454663.94696 
2454663.94788 



0.99986 
0.99970 
1.00028 
1.00031 
0.99885 



Note. — The full version of this table is 
available in the online edition, or by request 
to the authors. 



TABLE 2 
System Parameters of WASP-10 



Parameter 



Value 



68.3% Confidence Interval Comment 



Transit Parameters 
Mid-transit time, T c [HJD] 
Orbital Period, P [days] 
Planet-to-star radius ratio, Rp/R* 
Planet-star area ratio, (Rp/i?*) 2 
Scaled semimajor axis, a/R* 
Orbit inclination, i [deg] 
Transit impact parameter, b 
Transit duration [hr] 
Transit ingress or egress duration [hr] 

Other Orbital Parameters 
Semimajor axis [AU] 
e cos w 
esintt) 

Velocity scmiamplitudc K* [m s -1 ] 

Stellar Parameters 
M« [Mq] 

R* [R@] 
p* [pq] 

log [cgs] a 
[M/H] 

T ctt [K] 

Planetary Parameters 
M P [M Ju „] 

Rp [Rjup] 

Mean density, p P [pj up ] 
loggp [cgs] 

Equlibrium Temperature T c g(R, /a) 1 / 2 



2454664.030913 
3.0927616 
0.15918 
0.02525 
11.65 
88.49 
0.299 
2.2271 
0.3306 



0.03781 
-0.0453 
0.0228 
533.1 



0.75 
0.698 
3.099 
4.627 
0.03 
4675 



3.15 
1.080 
3.11 
3.828 
1370 



±0.000082 
±0.0000112 
-0.00115, +0.00050 
-0.00028, +0.00024 
-0.13, +0.09 
-0.17, +0.22 
-0.043, +0.029 
-0.0068, +0.0078 
-0.0075, +0.0098 



-0.00047, +0.00067 
±0.02 
±0.03 
±7.5 



-0.028, +0.040 
±0.012 
±0.088 
-0.0093, +0.0101 
±0.2 
±100 



-0.11, +0.13 
±0.020 
±0.20 
±0.012 
±50 



A 
D 
A 
A 
A 
A 
A 
A 
A 



B 

C 
C 

c 



B 
B 



A 
B 
U 
D 



B,C 

B 
B,C 

A 

D 



Note. — Note. — (A) Determined from the parametric fit to our light curve. (B) Based on group A 
parameters supplemented by the Y 2 stellar evolutionary models. (C) Based on our analysis of the C08 
RV measurements. (D) Reproduced from C08. 

a The logg* in the table is the value implied by the Y 2 stellar evolution models, given the measured 
values of T e g , [M/H], and p*. It is in near agreement with the value 4.40 ± 0.20 reported by C08, based 
on the widths of pressure-sensitive absorption lines in the stellar spectrum. 



